arXiv:chao-dyn/9502024vl 2 Mar 1995 


Chaos induced by Pauli blocking 


S. Drozdz 1 , J. Okolowicz 1 , M. Ploszajczak 2 , E. Caurier 3 and T. Srokowski 2 


1 Institute of Nuclear Physics, PL - 31-342 Krakow, Poland 
2 GANIL, BP 5021, F-14021 Caen Cedex, France 
3 Centre de Recherches Nucleaires, F-67037 Strasbourg Cedex, France 


ABSTRACT 

Dynamics of classical scattering in the system of fermions is studied. The 
model is based on the coherent state representation and the equations of motion 
for fermions are derived from the time-dependent variational principle, ft is found 
that the antisymmetrization due to the Pauli exclusion principle, may lead to hyper¬ 
bolic chaotic scattering even in the absence of interaction between particles. At low 
bombarding energies, the same effect leads to the screening of the hard, short-ranged 
component in the two particle interaction and thus regularizes the dynamics. 
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One of the central issues in the theory of complex dynamical systems is the 
classical-quantum correspondence for classically chaotic motion. Classically, chaos 
is a well defined concept but a way it manifests itself on the quantum level still re¬ 
mains a matter of controversy. One of the reasons of this controversy is the quantum 
mechanical symmetry related to the identity of particles. Of particular interest and 
importance in this context are the effects resulting from antisymmetry of an under¬ 
lying wave function for fermionic systems. Most of the physical objects, like atoms 
or atomic nuclei, providing empirical input for the above studies, are collections of 
fermions. The fundamental question then is in which sense does a fermionic nature 
of the basic constituents influence the structure of the corresponding classical phase 
space. And, in particular, does it increase or reduce the amount of instabilities 
leading to chaotic behavior? 

The model which seems best suited to incorporate the elements needed and thus 
for addressing such questions is the one based on the coherent state formalism. This 
formalism allows to define a classical limit JT] and, at the same time, to accommo¬ 
date the effects of antisymmetrization. Consistently, the i-th single particle state is 
described by the Gaussian wave packet : 



( 1 ) 


where Zi is the location of the center of gravity of the packet. These single par¬ 
ticle states are not mutually orthogonal and, thus, their overlaps n %] = {(j>i\(j)j) = 


exp [—(Zi* — Zj) 2 /(46 2 )] do not vanish. The A-fermion system is then described 


by a Slater determinant &(Z) = (A!) -1 / 2 det[0;( ?",,■)] whose norm N(Z) is equal to 
det[rijj]. In this way the state of the system is entirely specified by the parameters 
Z = {Z{ : i — 1,..., A} which for a proper treatment of the dynamics have to 
be considered as complex (Z = R + iP) J2J. The units are specified by setting 
h = b = m = Uuj = 1. The time development of the dynamical variables can be 
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determined by the time-dependent variational principle: 


dt($(Z)\(ij t - H)\<S>(Z))N(Z )-■ = 0, 


( 2 ) 


The resulting equations of motion take the form: 


* SiajpZjp 
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dH 

dZL 


( 3 ) 


where a,/3 — x,y,z. H is the expectation value of the many body Hamiltonian 
H(Z*, Z) = ($(Z)|Hj<f>(Z))iV(Z*, Z)~ l and the Hermitian matrix: 

d 2 
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is positive definite. Such a scheme constitutes a formal basis for various molecular 
dynamics approaches § B § and the Gaussian form of the wave packet proves 
appropriate for semiclassical studies of the structure of the phase space 0. 

In general, due to the antisymmetrization, the off diagonal terms in S ta jp do 
not disappear. Therefore, neither Z,- and Z* nor their real Ri and imaginary P, 
parts form the canonically conjugate variables. For the two particle system, an 
exact transformation can be performed to the new, canonically conjugate variables 
Wi (i = 1,2) These Wi cannot get closer than \/2 independently of the dif¬ 
ference between Z\ and Z 2 . Thus a topological hole, existing also for A > 3, 
corresponds to the Pauli forbidden region m- However, an explicit expression for 
the canonical variables is not known for A > 3 and, therefore, the equations (Q) will 
be solved in original variables Z. 

Schematic models 0,0,0. based on the scattering of a particle on the target 
composed of three particles at the corners of equilateral triangle located in the 
reaction plane, proved very instructive in studying various aspects of the collision 
processes. In this case, the scattering problem reduces to two degrees of freedom. At 
the present exploratory stage we perform an analogous study. The target is given 


3 




by a three fermion (Gaussian wave packets ([[])) configuration forming equilateral 
triangle of side equal to 4. We begin by entirely discarding the interaction term 
in the Hamiltonian in order to elucidate on the role of antisymmetrization itself. 
The Pauli forbidden regions manifest themselves by the strong increase of H(Z*, Z ) 
generated by the kinetic energy operator. Equivalently, the form of the energy 
surfaces in R depends on the momentum like variables P, as demonstrated in Fig. 1. 
Here the R dependence of the energy, as ‘seen’ by the fourth fermion located in the 
reaction plane (x,y), is plotted at two fixed values of P. With increasing P, the 
hilly structures in this effective Pauli potential become relatively smaller and the 
dynamics approach a free motion. Notice that a P dependence of the Pauli potential 
is an analog of the nonlocality for ordinary potentials. 

In a real dynamical process P changes in time and so does H . Nevertheless, 
the appearance of the three center structures suggests SB p~lH that for certain P 
values one may expect a strong sensitivity on the initial conditions. That this really 
happens is documented in Fig. 2 which shows, for the three different energies, the 
impact parameter dependence of the deflection angle 9 of the particle scattered off 
the target. The motion is initialized by setting the appropriate initial value of P x 
with P y = 0 and the impact parameter R y . The presentation of the results is based 
on the asymptotic values of the variables in the region where antisymmetrization 
is no longer effective and, thus, R , P are canonical conjugate. For energies either 
small (E = 0.2) or large (E = 0.6), compared to the height of the effective Pauli 
potential V m ^ (P = 0) = rna xH(R,P = 0) ~ 1/2 (see Fig. la), 9 depends on the 
impact parameter but this dependence is essentially continuous. At E = 0.4 and for 
R y between 0 and 1, one observes behavior characteristic of the chaotic scattering. 

Very interestingly, this scattering process carries all characteristics of the hyper¬ 
bolic chaotic scattering B. For this type of scattering, theory predicts 
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that the survival probability, i.e. number N(t) of trajectories remaining in the in¬ 
teraction region up to time t, follows: 

N(t) — N 0 exp[—A (1 — D)t\, (5) 


where A is the Lyapunov exponent and D is the fractal dimension. The Lyapunov 
exponent A, calculated from the growth rate of separation ratio 5/5q between the 


two neighboring (and long enough) trajectories [L|], is independent on what pair of 
scattering trajectories is used and, thus, A = 0.22 is well determined (Fig. 3a). Fur¬ 
thermore, the set of singularities seen in Fig. 2 for E — 0.4 possesses a well defined 
fractal dimension D. This one can conclude from calculation based on the uncer¬ 
tainty exponent technique |L5] which gives D = 0.591 (Fig. 3b). Finally, following 
the concepts of the transport theory [|T6|, by uniform random sampling of the whole 
interval of impact parameters one determines the survival probability N(t). In our 
case one observes (Fig. 3c) asymptotically exact exponential dependence, character¬ 
istic of the hyperbolic chaotic scattering, perfectly described by the above values of 
A and D. 

These investigations provide convincing evidence that the chaotic scattering we 
deal with is of the hyperbolic type where the set of singularities is connected to 
the existence of only unstable periodic orbits. In the present case, however, these 
structures take place for a system of noninteracting fermions and are entirely due 
to the antisymmetrization. Independently, it is interesting to notice the appearance 
of such uniform fractal set of singularities for a nonlocal problem. 

The most important result of the above analysis is that for a sufficiently dense 
system, the correlations resulting from Pauli blocking may lead to a significant mod¬ 
ification of the particle dynamics. In certain situations, they may even convert free 
motion in the gas of non-interacting particles into the strongly chaotic one. In gen¬ 
eral, this is more likely to occur when the mean kinetic energy of particles in the gas 
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is comparable with the height of the effective Pauli potential. At very low energies, 
the range of Pauli blocking for sufficiently dense system extends to such a size that 
different topological holes start overlapping and the dynamics becomes regular again 
(see upper panel of Fig. 2). In more realistic case of interacting fermions, this effect 
may thus screen out the short range components of the two-body interaction. This 
dynamical screening may be especially important for the hard, short-range interac¬ 
tions like the one between the nucleons. Classically, such an interaction generates 
strongly chaotic behavior. However, its appearance inside the Pauli forbidden region 
may restructure and even eliminate the corresponding irregularities. 

For a somewhat more quantitative illustration of this point, we present in Fig. 4 
the energy surface plots in the similar configuration as in Fig. 1. This time, how¬ 
ever, the constituents interact with the spherically symmetric, repulsive two-body 
interaction V(r) = Vo exp[—(r/ro) 2 ], with Vo = 10 5 / 2 and r o = 10“ 1//2 . For these 
parameters the interaction is comparatively hard and short-ranged. Upper panel of 
Fig. 4 corresponds to the situation with no antisymmetrization included. In (a) the 
side of the triangle equals 3 and in (d) it equals 4. In both cases the three hill struc¬ 
ture shows up and, consequently, the scattering will be chaotic in the corresponding 
energy intervals. Including antisymmetrization changes the picture completely (see 
Fig. 4b). Here, not only the energy is reduced by almost two orders of magnitude but 
also the three hill structure, previously responsible for chaotic behavior, disappears. 
Increasing momentum, slowly recovers the original shape of the energy surfaces, as 
is shown in Fig. 4 (c) and (f) for P x = 2y/l0, but they still remain about one order 
of magnitude lower. 

Returning to our introductory question we thus conclude that the fermionic 
nature of particles may drastically change the structure of the corresponding classical 
phase space. As the two extreme possibilities we identify the chaotic behavior in 
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absence of any interaction and the regularization of motion for strongly interacting 
particles. The first of them is very intriguing in view, for instance, of the fact 
that the hyperbolic chaotic scattering is considered as a classical manifestation of 
Ericson fluctuations 0- Equally interesting from the physical point of view is the 
second effect. The dynamical screening of the short-range components of the two- 
body interaction, is consistent with the success of the mean field concept in the 
region close to the Fermi surface for such a strongly interacting system as an atomic 
nucleus. In this connection, the space nonlocalities are identified |H| as a crucial 
element enhancing the nucleon mean free path. Many related questions remain a 
subject for further investigations. 

This work was partly supported by KBN Grant No. 2 P302 157 04. Two of us 
(S.D. and J.O.) wish to thank GANIL for the hospitality extended to them, when 
the work has been completed. 
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Figure captions 


Fig. 1 Energy surface plots for noninteracting configuration generated by the three 
fermions located at the corners of equilateral triangle as seen by the fourth 
fermion at momentum P x = 0 (a) and P x = a/TO (b). Triangle side length 
equals 4 and P y = 0. 

Fig. 2 Deflection angle 6 as a function of the impact parameter for three different 
energies under geometrical conditions corresponding to Fig. 1. 

Fig. 3 (a) Separation ratio between the neighboring trajectories for the scattering 
process at E = 0.4. Dots denote the dynamically determined values and the 
straight solid line represents a fit whose slope corresponds to the Lyapunov 
exponent A = 0.22. 

(b) Dependence of the fraction /(e) of uncertain pairs of trajectories (differing 
in 6 by more than vr/2) as a function of the difference e in initial values of 
the impact parameters. According to the uncertainty exponent technique the 
corresponding fractal dimension is D = 0.591. 

(c) Survival probability expressed as a number of the scattering trajectories 
remaining in the interaction region up to time t. The straight line represents 
the theoretically determined dependence. 

Fig. 4 Energy surface plots for the three fermion configurations forming equilateral 
triangle of side equal to 3 (a, b, c) and 4 (d, e, f) respectively: (a) and (d) 
represent the potential energy without antisymmetrization for the two-body 
interaction defined in the text, (b) and (e) the total energy with antisym¬ 
metrization, (c) and (f) the total energy with P = 2y/l0 
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